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We extend the theory of weak gravitational lensing to cosmologies with generalized gravity, de- 
scribed in the Lagrangian by a generic function depending on the Ricci scalar and a non-minimal 
coupled scalar field. 

We work out the generalized Poisson equations relating the dynamics of the fluctuating components 
to the two gauge invariant scalar gravitational potentials, fixing the new contributions from the 
modified background expansion and fluctuations. 

We show how the lensing equation gets modified by the cosmic expansion as well as by the presence 
of the anisotropic stress, which is non-null at the linear level both in scalar-tensor gravity and in 
theories where the gravitational Lagrangian term features a non-minimal dependence on the Ricci 
scalar. Starting from the geodesic deviation, we derive the generalized expressions for the shear ten- 
sor and projected lensing potential, encoding the spacetime variation of the effective gravitational 
constant and isolating the contribution of the anisotropic stress, which introduces a correction due 
to the spatial correlation between the gravitational potentials. 

Finally, we work out the expressions of the lensing convergence power spectrum as well as the 
correlation between the lensing potential and the Integrated Sachs- Wolfe effect affecting Cosmic 
Microwave Background total intensity and polarization anisotropies. 

To illustrate phenomenologically the new effects, we work out approximate expressions for the quan- 
tities above in Extended Quintessence scenarios where the scalar field coupled to gravity plays the 
role of the dark energy. 



I. INTRODUCTION 

In the recent years we reached a remarkable convergence in the most important cosmological pa- 
rameters describing the cosmic content as well the statistics of perturbations. The Universe is nearly 
geometrically flat, with an expansion rate of about 70 km/sec/Mpc, and structures grown out of a 
primordial linear spectrum of nearly Gaussian and scale invariant perturbations in the distribution of 
the energy density. About 5% of the critical energy density is made of baryons, while the remaining 
dark part is supposed to interact at most weakly with the baryons themselves, since we observe it only 
through its gravitational effects. The dark component appears to be 30% pressureless, like in Cold 
Dark Matter (CDM) scenarios, dominating the gravitational potentials perturbations which host visible 
structures like galaxies or clusters. The remaining 70% should be in some sort of vacuum energy, with 
negative pressure acting as a repulsive gravity, and responsible for a late time cosmic acceleration era. 
The case for this "concordance cosmological model" is now quite robust, supported by several indepen- 
dent datasets: the distant type la supernovae (hereafter SNIa 0,0), the Cosmic Microwave Background 
(CMB) anisotropies (see Q and references therein), the Large Scale Structure (LSS, 0, 0|) and the 
Hubble Space Telescope (HST @). 

The picture is clearly far from being satisfactory: in particular, without a better insight into the na- 
ture of the dark cosmological component, we cannot claim to have a satisfactory physical understanding 
of cosmology. The simplest description of the vacuum energy responsible for cosmic acceleration is a 
purely geometric term in the Einstein equations, the Cosmological Constant. On the other hand, while 
the CDM has a well established support from the theories beyond the standard model of particle physics, 
a Cosmological Constant providing the 70% of the critical density today raises two problems. The first 
is the fine tuning required to fix the vacuum energy scale about 120 orders of magnitude less than the 
Planck energy density which is supposed to be the unification scale of all forces in the early Universe. 
The second is a coincidence issue, simply why among all the small non-zero values of the Cosmological 
Constant, the value was chosen to be comparable to the critical energy density today. These questions, 
still largely unsolved, could be answered only if the concept of Cosmological Constant is extended to 
a more general one, admitting a dynamics of the vacuum energy, known now as the dark energy (see 
0,0111 an d references therein). 
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The simplest generalization, already introduced well before the evidence for cosmic acceleration 

nana, 

is a scalar field, dynamical and fluctuating, with a background evolution slow enough to mimic a con- 
stant vacuum energy given by its potential, providing cosmic acceleration. As soon as the latter was 
discovered, a renewed interest in these models appeared immediately 0, In particular, it was 

demonstrated how the dynamics of this component, under suitable potential shapes inspired by super- 
symmetry and super-gravity theories (see |l4j | and |15| . respectively, and references therein), can possess 
attractors in the trajectory space, named tracking solutions, capable to reach the present dark energy 
density starting from a wide set of initial conditions in the very early universe, thus alleviating, at least 
classically, the problem of fine-tuning [161 Il7| . The scalar field playing the role of the dark energy was 
named Quintessence. Its coherent insertion am ong the other cosmological components allowed to con- 
strain it from the existing data 0, 0, |2(], 0, H3, 0i H3 > as wen as to investigate the relation of 
the dark energy with the other cosmological components: the explicit coupling with baryons is severely 
constraine d by observatio ns [261, while the possible coupling of the Quintessence with the Ricci scalar 
[S3, IH, IH, Hfl, HH H3, IS3, EE [H, E3| and the dark matter [U [H [U , as well as the 

phenomenology arising from generalized kinetic energy terms |45l l46t \at\ have been extensively studied. 
The generalization of cosmology that we consider here concerns the gravitational sector of the funda- 
mental Lagrangian, admitting a gen eral dependence on the Ricci and Brans-Dicke scalar fields. This 
subject is interesting per se (see |48l |49| and references therein), and is receiving more attention after 
the discovery of cosmic acceleration, with the attempt to interpret the evidence for dark energy as a 
manifestation of gravity; this scenario has been recently proved to have relevant consequences for what 
concerns the dark energy fine-tuning problem mentioned above [40| . 

The gravitational lensing effect in cosmology is gathering a great interest becoming one of the most 
promising tools to investigate cosmological structures (see [5(j and references therein), and commonly 
thought in terms of strong and weak regime. The first one concerns highly magnified sources, gcncrically 
through the generation of multiple images of a background object by a single lens with the typical size 
of a galaxy. 

The weak gravitational lensing, which is the subject of the present work, produces weak amplification, 
generically through the distortion of the pattern of background light; moreover, it is generated by a 
large set of scales and objects, ranging roughly from non-linear structures like galaxy clusters to the 
large scale distribution of matter, still in linear regime. The weak lensing shear was detected recently 
by independent groups with astonishing agreement [5ll l52l l53l l54l l55| . Although the precision of such 
measurements does not allow to constrain different cosmological models, the planned observations will 
become certainly a crucial tool to investigate the behavior of dark matter and energy during the structure 
formation process and at the onset of cosmic acceleration j5|J . 

In particular, several authors considered the possibility to investigate the dark energy component through 
the weak lensing, with different approaches ranging from shear distortion of background galaxies from 
clusters of gala xies to the weak lensing power causing a non-Gaussian pattern into the CMB anisotropics 
[571 l58l l59l l60l l6lL 162. l63l l64j . The reason of this interest is the timing: the structure formation, and 
the weak lensing carrying its physical information, occurs at an epoch which overlaps with the onset 
of cosmic acceleration; by virtue of this fact, it is reasonable to expect a good sensitivity of the weak 
lensing effect to the main dark energy properties such as the equation of state and its redshift behavior. 
These studies are entering in an higher level of sophistication: the first outcome of N-body simulations 
in several dark energy scenarios have been published recently [6^, |6(| [63 , mainly studying the impact 
of the background rate of expansion on the internal parameter of structures. The implementation of a 
light ray tracing technique through those structures allows to check numerically the weak lensing pattern 
produced by mildly and full non-linear density perturbations |68l ] . 

The theory of weak gravitational lensing in ordinary cosmologies, i.e. made by radiation, dark matter 
and cosmological constant, is known and well established, but a comprehensive treatment in cosmologies 
with generalized theories of gravity still lacks in literature. The aim of the present work is to fill this 
gap, by providing the community with the recipe to interpret the weak lensing observations in a more 
general context. We follow the harmonic approach to weak lensing |69|| and the treatment for generalized 
cosmological scenarios including cosmological perturbation |49| . already exploited for investigating the 
effects of the explicit coupling between dark energy and gravity (see [4(| and references therein). 

This work is organized as follows. In Section II we describe the cosmological models we deal with, for 
background and linear perturbations. In Section III we write the generalized Poisson equations relating 
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the fluctuations in metric and cosmological components. The derivation and discussion of the lensing 
equation and weak lensing potential are in Section IV and V, respectively. Finally, in Section VI we 
draw the conclusions. 



II. GENERALIZED COSMOLOGIES 



We shall consider a class of theories of gravity whose action is written in natural units as 



S = Id x^J—g 



(1) 



where g is the determinant of the background metric, R is the Ricci scalar, uj generalizes the kinetic term, 
and £fluid includes contributions from the matter and radiation cosmological components; k = 87rG» 
plays the role of the "bare" gravitational constant, [34[. Here as throughout the paper Greek indices 
run from to 3, Latin indices from 1 to 3. 

The usual gravity term R/l6nG has been generalized by the general function f /2k |3lL l49l|. Note that 
the formalism adopted is suitable for describing non-scalar-tensor gravity theories, i.e. without (j>, where 
however the dependence on R is generic. 

A. Einstein equations 

By defining F(<j),R) = (l/n)df /dR, the Einstein equations G M „ = kT^ v take the following form: 

where again one can recognize a part depending on the fluid variables, and a part relative to the non- 
minimally coupled scalar field of the theory plus a contribution arising from the generalized gravity 
coupling represented by the function 1/F; for practical purposes we will render this splitting explicit 
rewriting T^ u as 

T^ = yTlf A + T^ . (2) 

Note that in our scheme the generalized cosmology term is active also if gravity is the same as in general 
relativity, and a minimally coupled scalar field, like in the Quintessence models, represents the only new 
ingredient with respect to the ordinary case. This is true also in the limiting case where the scalar field 
reduces to a Cosmological Constant. Moreover, notice that is not conserved if gravity differs from 
general relativity; nonetheless, contracted Bianchi identities still hold and ensure 

(T;ry = o, ( 3 ) 

which allows to derive the equations of motion for the matter variables only, leading to a remarkable 
simplification 49] . The expression for the stress-energy tensor relative to the scalar field which is 
conserved in the generalized scenarios described by the action £Q) must include the term accounting for 
the interaction with the gravitational field |33]. It is also worth to notice how the equations (J2J) get 
simplified if the function / is a product of R times a function of the field only: 

f 

— = RF . (4) 

K 

In the following, we will refer to this class of cosmologies as Non-Minimally Coupled (NMC) models. 
We willingly keep this work as much general as possible, due to the large variety of scenarios covered by 
the actionQ] We shall only consider the Extended Quintessence (EQ |33Ll38l p scenario as an example to 
illustrate the new aspects of the weak lensing process with respect to ordinary cosmologies; in that cases, 
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the field <j), non-minimally coupled to gravity, also represents the dark energy, providing acceleration 
through its potential V. Specifically, the original works considered a NMC model defined as 

F(<P) = -+^\ (5) 

K 

and an inverse power law potential V(<f>) — M i+a /(jf 1 providing cosmic acceleration today. The con- 
straints from solar-system experiments force the correction to the gravitational constant to be small 
in this specific models. Therefore it is suitable to make approximations to illustrate a sort of first or- 
der variation of the weak lensing in generalized gravity theories with respect to the case of ordinary 
cosmology. 



B. Background 

We will write the unperturbed Friedmann Robertson Walker (FEW) metric in spherical coordinates 
as: 

ds 2 = a 2 (r) (- dr 2 + - ^ Rr2 dr 2 + r 2 dnj , (6) 

where K is the uniform spatial curvature of a spherically symmetric three-space, d£l is the metric of the 
two-sphere, and r stands for the conformal time variable, related to cosmic time by the usual relation 
dt = a(r) dr. 

The energy-momentum tensor can be recast in a perfect-fluid form: 

T^v = (p + pju^Uv + pg fiV ; (7) 

the corresponding background energy density and pressure are easily computed to be: 

1 ( w J/2 RF-f/n 3HF'\ 1 

P = j I Pfluid + ^ + 2 + ~ ) = F Pfiuid + Pec ' ( > 



1 ( ^ j.,2 RF-f/n y t F" HF'\ 1 
p= F\ Pa * d+ 2rf 4> 2 V + ^2+—)=F Piiuid+p z c > (9) 

the prime denotes differentiation with respect to conformal time and TL is the conformal Hubble factor 
a' I a. As above, p gc and p gc do not obey the conservation law in ordinary cosmologies, p' gc + 3Ti(p gc + 
Pec) ¥= 0. ^ ~ ... 

In FRW cosmologies the expansion equation reads 

H 2 =a 2 p~K, (10) 

and it cannot be solved directly due to the appearance of Ti. in p gc , which is explicit in the last term but 
is also contained in RF — f/n through 
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R=-(H + H 2 ). (11) 



Note that this is true also in theories where / = f(R) and no scalar field is present. 
On the other hand, NMC scenarios admit a formal solution, which is 



]jl (? ) 2 + 1 ( a2pfluid + 1^' 2 + a2y ) K ' ( 12 ) 

where we have selected the expansion solution with positive TL. Note that the dependence of the comoving 
distances r on the redshift z = 1/a— 1 gets also modified, according to (|12f> : 

o H{z) 
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We generically indicate the single components in the fluid with x. Since Tg uid is conserved, energy 
density, pressure, equation of state and sound velocity, defined as 

Px = ~T$ X , p x = 1/3 T\ x , w x = p x /p x , c 2 s = p'Jp' x , (14) 

give rise to conservation equations having the familiar form: 

p x + 3Hp x (l+w x )=0 . (15) 

The last ingredient is the Klein-Gordon equation for the evolution of the field, which is substantially 
different from the case of ordinary cosmologies: 

0" + 2U4>' + (u,^' 2 - a 2 ^ + 2a 2 V^j = . (16) 

As we stress in the next Section, the relevant changes with respect to the standard picture are represented 
by the change in time of the function /. In EQ scenarios the dynamics of the field possesses two distinct 
regimes. At low redshift, the behavior of the energy density coincides with the corresponding one in 
the tracking trajectories in ordinary quintessence models, linked to the potential exponent as — 
—2/(2 + a). At high redshift, generally much earlier than the epoch of structure formation, eventually 
the effective potential coming from the non-minimal interaction wi th g ravity takes over (R— boost), and 
imprints a behavior w<p = — 1/3 for the quadratic coupling (JSJ |U l40j ] . 



C. Linear cosmological perturbations 

We will describe the linear cosmological perturbations in the real as well as in the Fourier space. 
For this reason, we follow a notation close to the one introduced recently by Liddle and Lyth 70], 
which allows not to explicitate the Laplace operator eigenfunctions when working in the Fourier space, 
minimizing the formal changes needed to go from the real to the Fourier space and viceversa. See jz3 
for the usual formulation cast in the Fourier space. 

The general expression for the linear perturbation to the metric © can be written as |70| 

ds 2 = a 2 (r){-(l + 2 A)dr 2 - B^rdx 1 + [(1 + 2D)8 rj + 2E ij \dx i dx i }, (17) 

where the function Eij is chosen to be traceless in order to uniquely identify the non-diagonal spatial 
perturbation. 

It is usual to further decompose the above quantities Bi and E^ into pure scalars (S), scalar- type 
(S) and vector-type (V) components of vectors, scalar-type (£>), vector-type (V) and tensor-type (T) 
components of tensors, according to their behavior with respect to a spatial coordinate transformation: 

A= A s ; (18) 
B t = Bf + BY; (19) 
E ij = EZ + EYj+Efii (20) 



where 



Bf = V,B, Eft = ( ViVj - -%V 2 ) E , (21) 



1 

3' 

V*BY = 0, EY = \ {V t E 3 + VjEi) , E\ v = , (22) 
Ej T = 0, V*£g = 0, (23) 



and E, B are scalar functions. Our notation slightly differs from the one of Liddle and Lyth which 
is given in the Fourier space: the quantities E, Ei and B we use here correspond to the original ones 
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divided by k 2 — kik 1 and k, respectively. In the linear theory the different types of perturbations evolve 
independently of each other and can thus be treated separately. 

An analogous decomposition can be performed for the stress-energy tensor, whose expression up to the 
first perturbative order is: 

T°=-(p+<5p), (24) 
Tf= (p+p^Vi-Bi) , (25) 
n= -(P + PW , (26) 
fj= {p + 5p)5)+pW j , (27) 



where the fluid velocity Vi and the anisotropic stress 11^ can be split as 

Vi= vf + vY; (28) 

HJ - ng + nl + ng 



n ii = iig + nK+n^ (29) 



with the same properties of their metric counterparts: 



ViVj 



r%v 2 ) n , 



vtf = o , n$ = - (v^ + v^) , n} v = o 
vn, = vtf". = o , 



(30) 

(31) 
(32) 



where again II is a scalar quantity and the same differences of our notation with the one by Liddle and 
Lyth hold here for the stress-energy tensor perturbations. 

In this paper we will take into account only density (i.e. scalar-type) perturbations. The reason is 
that also in the generalized scenarios we consider here they play the dominant role. In fact, as we stress 
in detail in the following, the scalar-tensor coupling does generate a non-null anisotropic stress already 
at the linear level, but that is of scalar-type only, and therefore does not act as a source for gravitational 
waves. 

We work in the so called conformal Newtonian gauge, where the non-diagonal perturbations to the metric 
are set to zero: 

B = E = 0. (33) 

Furthermore, we will rename the lapse function A and the spatial diagonal perturbation D after the 
widely used gauge- invariant potentials 

A -> * ; D -► $ . (34) 

The line element used through the rest of the paper will therefore be 

ds 2 =a 2 [-(l + 2*)dT 2 + (l + 2$)^ 2 ], (35) 

where dl 2 is the unperturbed spatial length element from ©. 

We now write down the main equations driving the evolution of the perturbed quantities defined above. 
For each fluid component, the evolution of the scalar perturbed quantities can be followed through the 
dynamical variables S x = Sp x / p x , v x , 5p x , W x , defined in terms of the stress energy tensor as 

S Px = ST$ X , 6 Px = 1/3 8Tl x , V l v x = -5T*J(p x + Px ) , ^V.V^TI, = 5T^ 1 . (36) 
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Note that from now on we do drop the subscript S meaning that we always treat scalar cosmological 
perturbations, unless otherwise specified. In the Fourier space, the equations for S x and v x take the form 

S' x - 3Hw x 5 X = k 2 (l + w x )v x - 3(1 + w x )& - 3Hw x p x 8p x ; (37) 



< + H(l - + = - j^j - * + | (l - f ) JL.II. , (38) 

while Sp x and depend on the particular species considered. The perturbed Klein-Gordon equation can 
be written in terms of two equations formally equivalent to (|37I38I) by building the conserved expression 
for the perturbed energy density, pressure and anisotropic stress perturbations |37|. Their combination 
leads to the Klein-Gordon equation at first perturbative order: 



1 /UJ 



uj / ,<f) 



(¥' - 3$'W + f + W + —<P' 2 ] * + ^ 
\ uj J 2a. 



-a 2 f^/n + 2a 2 V 4 



2u> 



:ujk i 



dR 



SR. 



(39) 



III. GENERALIZED POISSON EQUATIONS 



In this Section we work out the equations relating the stress-energy tensor perturbations to the scalar 
metric gravitational potentials. The latter, together with the background cosmic geometry and expan- 
sion, determine entirely the lensing process. 

We start writing the generalized expression of the density fluctuation ^| : 



Sp = p ■ 8 



ST, 



H5F' 
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1_ 

F 



Sp 



'fluid 



U) 

— ' 

■ft A 



u>,4> 



f, 



2V e , 



p + 3p 



1 



SF 



&F' 



(40) 



We can focus on two main aspects of the generalized expression above, playing the major role into the 
generalization of the Poisson equation: the contribution from the field fluctuations 8(f> and the 1/F term 
in front of the expression for Sp, which acts as an effective time varying gravitational constant. As we 
shall see below, the latter is the relevant effect in typical Extended Quintessence models. 
The Poisson equation relates the fluctuations in the time-time component of the metric to the usual 
combination of density and scalar-type velocity perturbations, named A, whose expression in Newtonian 
gauge is |7l| : 



A = 6 -3Hwv . 



(41) 



We follow as much as possible the notation of earlier works [Hij- The SGq = ST® equation can be cast 
in such a way that formally it coincides with the case of ordinary cosmologies. In the Fourier space it is 



3 A') <I> = 3 A j H 2 + — 



(42) 



so that we can exploit our distinction between fluid and generalized cosmology terms. Note that the 
Hubble expansion rate is evaluated with respect to the ordinary time, H = a/a = a 1 /a 2 = Ti/a. By 
using equations IjHIlUfl . we can write H 2 



-^fluid 



H 2 C , where 



H 2 



fluid 



1 Pflui 

F 3 



F 



n 0m (i + zf + n 0r (i + zf + (i - no) (1 + zf 



(43) 
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and f2omj ^Or and f^o are the contribution to the present expansion rate from the matter, radiation and 
the total density respectively, while Fq is the actual value of the gravitational coupling strength, and can 
be replaced with 1 /8ttG. It is important to note that H^ uid is linked to the energy density of the fluid 
components, all them but the scalar field, but it contains a most important generalization, represented 
by the Fq/F term, which plays the role of a time dependent gravitational constant into the Friedmann 
equation. Moreover, using the relation K/a 2 — H 2 (£Iq — 1)(1 + z) 2 , we included the effect of the spatial 
curvature into -f/fluid- The expressions for H 2 C and fl gc can be easily obtained by making use of the 
equations (|8I1U|) : 

H* C =P SC , «gc = ^. (44) 

Starting from (|42|l . let us write down the relation between the power spectra of <&, P$ = k 3 & 2 /2n 2 , and 
the one of A: 

P 9 = S.(**i) i (i 3(i-no)zr aN 1 



4 V k J V k 



2 



^o 0m (i + z) + ^n 0r (i + zf + (f - n„) + (TT^ gc 



-| 2 

Pa • (45) 



This is the equation which generalizes the link between time-time metric fluctuations with density and 
scalar-type velocity perturbations. The new effects arise from the scalar field contribution, encoded in 
H 2 C , the 1/F term behaving as a time dependent gravitational constant, as well as from the fluctuations 
of the scalar field, contained in A, both in bp and v. 

Let us check the most relevant corrections to the quantities above in EQ models we take as reference 
I13,|3^|, where the scalar field fluctuations Scj) play a minor role [13, and the overall geometry is assumed 
to be flat, K = 0. In these condition, the most important correction is represented by the 1/F term, 
effectively representing the time variation of the gravitational constant. The expression for F in can 
be conveniently rewritten as 

F= 8^G + ^ 2 ~^ } ' (46) 

to make explicit that at the present F = 1/ 8nG. The observational constraints 0] usually are 
expressed as bounds on the quantities 



1 dG ^ 1 dF _ 2 d<t> F 
G~dt ~ ~F~dt ~ 0~dt ' UJJBD ~~F 2 4£ 



WJBD — "^2 = 77T72 ' (47) 



calculated at the present time, where the ~ sign above is due to the slight difference with the gravita- 
tional constant measured in Cavendish like experiments 1 3 Jj j and lojbd is the usual Jordan-Brans-Dicke 
parameter, which usually implies the strongest constraint. Typically [j^HH the correction to the I/87rG 
term is small, so that 

j ~ 8nG[l - 8ttG£(0 2 - <&)] . (48) 

Moreover, in tracking trajectories with inverse power law potentials, the field approaches <f>o of the order 
of 1/VG from below, being generally much smaller than that in the past, when 1/F freezes to the value 
87rG(l + 87tG£0q). The magnitude of the correction is therefore 



SttG^ 2 = = ^ J — ■ (49) 

^OJJBD V UJBD 

note that as ujjbd approaches infinity, recovering general relativity, the correction may still be relevant 
depending on the values of £ or </>q. The gravitational potential receives a contribution which is 

5<S> = -8ttG£(0 2 - ^)$ fl uid+0 ■ (50) 
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The subscript ()fl u id+0 represents all the terms coming from the fluid quantities as well as the scalar 
field ones from the minimal coupling, i.e. not involving F explicitly: 



k $fiuid+0 = 47rGA flui d+0 J ffo Pom (1 + z) + fl 0r (1 + zf + fy] 



(51) 



Note also that if the trajectory is tracking, with an almost constant equation of state iu^, the expression 
above reduces to 



fc 2 $fluid+0 = 47rGA fluid+0 H o 2 [r! Om (1 + zf + f> . r (1 + zf + fi 0o (1 + zf (1+w ^} ; 



(52) 



in order to keep the notation simple, we drop such subscript in the following, always meaning that it 
is there when discussing approximate expressions in EQ models. Similarly, the gravitational potential 
power spectrum gets an extra contribution which at the first order in the correction to 1 /8nG is 



6P® = -16ttG£(0 2 - <pl)P$ 



where 



Ho 

k 



^Orn (1 + Z) + Qor (1 + zf + Sl^ (1 + z) 



Pa 



(53) 



(54) 



In general, another important effect which arises in cosmology from the generalization of the underlying 
gravity physics is represented by the relation between \? and The difference between / and R, which 
may arise from a scalar-tensor coupling as well as a non-standard dependence of / from R itself, gives 
origin to tidal forces exciting the anisotropic stress of scalar origin 49] ; in the Fourier space its simple 
form is 



k 2 6F 



(55) 



and implies a shift between the two gauge independent scalar metric perturbations, which in our gauge 
takes the form 



a 2 SF 

* + * = -ppn = --, 



(56) 



where the last equality holds if the anisotropic stress is due to the generalization of gravity and is 
does not come from matter or radiation. This is an important new aspect of generalized cosmologies 
which implies a change in almost all the equations describing the weak lensing effect in cosmology, to be 
discussed next. For this reason, it is convenient to give a name the "J + <& combination, valid both for 
the real and the Fourier space: 



S = * + $ . 



(57) 



One has S = in ordinary cosmology, and 3 = —5F/F in the generalized scenarios of interest here. As 
we already stressed, S is excited both by a scalar-tensor coupling like in EQ models, and a generalized 
dependence of the gravitational Lagrangian term on R] its expression in terms of / is 



dR 



d 2 f , 
d(f>dR l 



dfi 

dR 2 



SR 



(58) 



While in the first case the correction is small because of the smallness of the scalar field fluctuations 
[3ll l37l] . the contribution from the second term has not been investigated yet. 



IV. LENSING EQUATION 



Before working out the generalized expression for the weak lensing power spectrum it is necessary 
to reconsider the lensing equation to track the new effects coming from background dynamics and 
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perturbations. We follow the approach by [7q. deriving the photon trajectories as solutions with ds 2 = 
of the geodesic equation for the metric i|17|) ; the lensing deflection around a given direction in the sky is 
described introducing new angular coordinates 9 X and 9 y , defined as 

9, x = 9 cos ip; 9 y = 9 sin ip , (59) 



where 9 — ^J9 2 + 9 2 and (p are the polar coordinates in the (9 x ,9 y ) plane. 
Furthermore, we perform a change of the radial coordinate: 



X(r) = — = arcsin \f~Kr (60) 
v K 



such that the spatial background metric takes the more readable form: 

d f = dx 2 + sin 2 VZ X (dg 2 + d92) (61) 
K y 

This notation is convenient because the weak lensing hypothesis immediately reflects in the condition 

9 <C 1 , (62) 

which allows to write the geodesic equation at first order in the deflection angle, in addition to the usual 
linear approximation for metric perturbations. The geodesic equation is indeed 

d 2 r a aP f 1 \ dr* dr" 

^ = -^^-^J^, (63) 

which for a — and a = 1 gives dr = l/a 2 d\ and dr = l/a 2 d\; substituting these expressions into the 
perturbed equation for the angular part we get 

d 2 9 x K „ , . / — cos \J~Ky d9 x , . 

-rf = 5-7^-^-(*- $ -2\/^ — 64 

dr 2 sin 2 V^X sin Vif X rfr 

for the angular coordinate and another one for (9 y , formally equivalent with x —* H- 
In terms of S, defined in the geodesic equation assumes the familiar form [50,173 plus the pertur- 
bation coming from the anisotropic stress (i stands now for x or y): 



dr 2 s in 2 \f~Kx * sin \[K\ dl ~ 



(65) 



The effects coming from the modified cosmological expansion are encoded in \, through the modified 
dependence of the distances r given by (|12I13|) with respect to the redshift z. 

In terms of the comoving displacement from the polar axis, Xi = 9i sin V^x/V^, the equation (|65|) 
simply reads 

x'l + Kxi = y^~^ =2di$-diZ, (66) 
axi 

where the first term on the left hand side describes the tendency of two nearby rays to converge, diverge 
or remain parallel according to the geometry of the universe, while the right hand side accounts for the 
lensing effect due to the metric perturbations. The general solution of this equation is 



x . + A SmV * X +Bcos^Kx - - f X d X 'd^(n, X ') ~ $(n,x')] 
V K Jo 



K 



X dx'd i [2$(n >X ')-E(n,x')} Sm ^ x ' , (67) 



K 

where A and B are integration constants, and the position x on the light cone is completely specified by 
the line of sight direction n and the generalized radial coordinate x- 
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V. WEAK LENSING 



It is convenient to begin with the comoving separation between two lensed rays, starting from the 
same point, one in the direction of the polar axis and the other one in a direction n, on a source plane 
at distance Xs', in terms of the angular separation Qi in the direction i one has 



Xj (Xs) = (Sij - ipy ) Sm ^ X 9 l , (68) 



where ipij is the distortion tensor: 



1 f x ^ , sin VKx' sin V~K( x ~ x') 
VK Jo sin \f~Kx 



ipij{n,X) = —fT7 I d X — Oi9,-[*(n,x ) - )] 



1 f x ,sinV^x'sinVX(x- x') a « r '\ , -c- 'm 

= ~ri7 d * : — 7^ a i d i [-2*(n,x)+a(n,x)]- (69) 

V A Jo sm y Kx 

Note that at first order the effect of cosmological perturbations can be computed along the unperturbed 
trajectories, which corresponds to neglect the difference in the relative deviation of two lensed and 
unlcnsed rays inside the above integral. 

The components of ipij are usually interpreted in terms of the shear 7 = 71 + 172 and of the effective 
convergence «, respectively identified as 

7i = ^(^n _ ^22) ; 72=^12; 

«=|(Vii+V^)- (70) 

Moreover, as we usually deal with lensing phenomena from a multiplicity of sources, the distortion tensor, 
and thus the projected potential, are usually meant to be integrated over the possible source distances: 

= J dxg(x)i } ij(n,x), (71) 

where g(x) is a normalized function describing the distribution of the relevant sources. When one 
considers the effect of lensing on the CMB, the source distribution may be replaced by a delta function 
at the last scattering surface. 

Let us evaluate the correction to the distortion tensor in EQ models. The contribution from S is negligible 
in these scenarios, since it arises from the scalar field fluctuations 5<fi, yielding a correction which is small 
with respect to the one coming from the variation of 1/F 37]. Therefore, in flat cosmologies and at a 
given x, the correction to ipij(n,x) is due only to the shift in the gravitational potential, represented by 
JSOJ: 

5^j(n, X ) = 167rG£$ /* d x ' X ' ix - X ° (€ - l) d t d^(n, x '). (72) 
Jo X \<Pq J 

Note however that integrating in the x variable, although convenient in order to minimize the formal 
corrections to ipij, hides the effect of the varying gravitational constant on x itself; in flat cosmologies 
the latter coincides with r given by l|13l) . which has to be corrected as 

mi*- 1 )- (73) 



as it can be easily verified since H oc 1/ VF. 
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A. Generalized lensing potential 

The lensing equation is often rewritten in terms of the projected potential </>, defined through the 
relation 

1>H = T . % d.^- (74) 
(sin V-rvx) 



Indicating the radial coordinate distance with D{x) = siny K\/y/ K, we get 

4>(n) = jT~ d X D( X ) x) - $ (n, x)] ^ ^^^ (x') 



° dx^(x) [-2*(ft, x) + S(n, x)] /*" ^^pp5(x / ) , (75) 

where Xoo stands for the comoving distance at infinite redshift. By defining the integral involving the 
source distribution as 

g'(x) = D( X ) dx' ^'^ gix') , (76) 

the lensing potential l|75(l takes the compact form 

<f>{ti) = [ °° d X g'(x)[y(n,x)-Hn,x)}= f °° d X g'{x) [~2*(fhX) + 3(ft,x)] • (77) 
Jo Jo 

The expression above acquires several new contributions in the generalized scenarios of interest here. The 
modified background expansion affects the angular diameter distance as well as the effective gravitational 
constant; the perturbations get new contributions from the field fluctuations affecting $ and exciting 
the metric fluctuation mode represented by S. 

In EQ models, if the integration is made on the variable X , the main correction is due to the time 
variation of the effective gravitational constant: 

ty(n) = 167rG£$ £°° d X g'{x) (78) 

In this equation, and in Eqs. Q86 fl .H9i p .H94 p . the quantity $ is meant to be $fl u id+0 as defined in Eq. 
(EDI . 

We need now to track these effects into the angular power spectrum of the projected lensing potential. 
That is defined as usual as 

= ( \(t>i m \ 2 ) , <km = -2 / rffift <f>(n) Y lm (n) (79) 

where the —2 is purely conventional in order to keep the notation consistent with earlier works |69| . 
One needs now to expand the metric fluctuations in the Fourier space with respect to the position 
x = r ■ n. The expansion functions are just the eigenfunctions YAx) of the Laplace operator in curved 
spacetime, defined in general in curved FRW geometry |7lj . Their radial and angular dependences are 
further expanded in ultra-spherical Bessel functions ui and scalar spherical harmonics, by exploiting the 
relation 

Y^)=ATrY J i l n{kx)Y lm {k)Yi m {x), (80) 
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where k and x denote the modulus of the corresponding vectors. By using the completeness of the 
spherical harmonics, and the fact that x coincides with the radial distance -D( X ), the final expression for 

4>lm IS 



dxg'(x) / d 3 k [2$(fc, X ) - H(fc,x)l i l ui(kD(x))Y lm (k)Yi m (n) 



(81) 



The lensing potential angular power spectrum Q79J1 is therefore 



32 



dxg'ix) / dx'g'ix') / d 3 k / d 3 fc' uz [fcD( x )] Ui [A' J D( x ')]^(fc)y^(fc') 



(<s>(k,xMk',x'y - ls(k,xMk',x'T - ~s(#,x7*(fc.x) + Js(£, x )s(fc', x 'r> 



(82) 



Assuming that the statistical average above eliminates the correlation between different Fourier modes, 
as well as the dependence on the direction of the wavenumbers, 

(A(k, x ) B(k', x')) = (A(k, x )B(k, x '))S(k - %') , (83) 
where A and B represent either ^ or S and they are meant to be ensemble averaged, one finally gets 



32 

7T 



dxg'ix) / dx'g'ix') / k'dkmikDix^mikDix')} 



*)*(*, x')) - (E(k,xWk, X ')) + i (S(fc, X )S(fc, X ')) 



(84) 



It is also useful to write down explicitly the equivalent form of the expression above which contains the 
gravitational potentials only: 



£,4>4> _ 32 



dxg'ix) / dx'g'ix') / k'dkmikDix^mikDix')} 



^(k,xMk,x')) + l($ik,xMk, x ')) - i(*(fc, x )$(fc, x ')) 



(85) 



putting in evidence the correlation between ^> and $. 

From this expression we can easily infer the main correction to the lensing potential angular power 
spectrum arising in EQ cosmologies, using l|53|l : 



512GC0O 



dxg'ix) / dx'g'ix') 



1 



o 



fc 2 dfc Uz [fcD( x )] Ui [fci)( x ')]($(fc, X )$(fc, X ')). 



(86) 



Note that the numbers here conspire to yield a quite large factor in front of this expression, which may 
render the correction above relevant even for values of the product G£</>§ as small as 10~ 3 . 

In the following we further specialize our results computing the generalized expression of some quantity 
particularly relevant for observations, as well as their main corrections in EQ models. 



B. Convergence power spectrum 

The convergence, represented by the trace of the distortion tensor, is usually used as a main magnitude 
of the weak lensing distortion. Indeed, in the weak lensing hypotheses, it coincides with the shear 
power spectrum (see |50|) which has been recently observed from the distortion induced in the shape of 
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background galaxies in the optical band 0, |5^, EH Ell • 

The expression we need to compute is given by 



Keff 



-(^11 +-022) 



(87) 



which shall be averaged over the source distribution as usual. We get 

Scff = \J^° d x 'g(x') f°° D{X) d{x') ^ * dl [ * {fl ' x) ~ $( "' x)] = 

1 f Xa ° f x ™ D(y) D(y' - y) 

= 2 1 dx ' 9ix ' ] J x D§) - y0i[- 2 g(".x) + B(n,x)] 



The two-dimensional Laplacian appearing in this equation can be safely replaced with its three- 
dimensional analogue (see |50( and [76J for a numerical check of this point). Once this substitution 
has been made, we can expand the generalized gravitational potential in Fourier harmonics transform- 
ing with respect to the spatial point h ■ Xj an d transform the Laplacian in a multiplication by (—k 2 ) : 



1 



Kcff 



2(2tt)3/2 J 



dxg'(x) / d 3 kk 2 [2$(£,x)-S(fc, X )]Yg(a?). 



(89) 



Comparing this expression with the ones for the lensing potential power spectrum 1)84185(1 we can imme- 
diately infer the result: 

Pk(l) = - J dxg'(x) J d X 'g'(x') J dkk s Ul (kD(x))ui(kD(x'))- 

mk,xMk, x ')) - (S(fc,x)$(fc,x')> + \(Z(k, X )Z(k,x')) 

X ~ dxg'(x) I*™ dx'g'ix') I dkk s ui(kD(x))ui(kD( x '))' 
o Jo J 

i(*(fc,x)*(fc,x')) + \{<s>(k,xMk,x')) - l{*(k,xMk, x ')) 



(90) 



In EQ cosmologies, using again (|53|> . one finds 



5P h {l) = -USG^l 



dxg'(x) dx'g'ix') l^-i 



dk k 8 Ul (k D( X ))ui(k D( X '))mk, xWk, x')> , 



(91) 



where the correction to x must be taken into account following l|73[) if the integration is made on the 
redshift. 



C. Lensing in the CMB signal 

The lensing potential correlates significantly with secondary anisotropics of the CMB, because it arises 
at the epoch of structure formation; here we generalize the lensing cross-correlation with the Integrated 
Sachs- Wolfe effect (ISW, sec 69] for a comparison with the case of ordinary cosmologies). 
The latter can be represented in terms of temperature fluctuations as 



e 



ISW 



(n) 



d X [Hn,x)-^(n,x)} = - / d X [2$(n, x) - S(n, x)] • 



(92) 
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Note that, in order to avoid confusion with the integration variable in this paragraph only we will 
denote with an overdot the derivative with respect to conformal time. 

Again using the expressions for the lensing potential power spectra (|84I85|) . and making use of the 
statistical independence of different Fourier modes, we are able to write immediately the cross-correlated 
spectrum: 



Q0<t> _ ° 



^ ' dxg'ix) I*™ d X ' I k 2 dk Ul (kD(x)) Ul (kD( x ')) 



($(fc, x )$(fc, x ')> - \m,x)HKx')) - \m,x'mk,x)) + \m, X )Hk,x')) 

2 /X ~ dxg'ix) f X °°dx' I PdkmikDixVmikDix'))- 



(*(fc, x )*(fc,x')> + <<f(fc,x)$(fc,x')> - <*(fc, X )$(fc,x')> - (*(k,x'Mk, x )) 



(93) 



The main correction in EQ cosmologies is obtained again by using (|53(1 . Interestingly, the time derivative 
reintroduces a term proportional to |$| 2 : 



8C^ = -128G^ 2 /*~ d X g'(x) /*~ d x ' (t 1 
Jo Jo Wo 



k 2 dk mik D( X ))u,(k D(x'))($(k, xMk, x'))- 



256G£0g / dxg'ix) / d x '- 



k 2 dk mik D( X )hi(k Dix'))(Hk, xMk, x')) ■ (94) 



These expressions can be further simplified noticing that for lensing on the CMB signal the source 
distribution is well represented by a delta function at the last scattering (LS) surface; thus the averaging 
function g'ix) can be written as 



't \ r>t \ f X °° a ' D (X' ~X) , , Dix)D(XLS-x) 

gix) = D ix )J x d x ^ r six-x,s)= d(xls) 



(95) 



VI. CONCLUSIONS 



The weak lensing in cosmology is one of the most important tools to investigate the structure formation 
process, and in particular the mechanics of the dark cosmological component, which represents almost the 
95% of the cosmic budget according to the recent measurements. Since the onset of cosmic acceleration 
occurs at the same epoch of structure formation, the weak lensing is most promising in order to gain 
insight into the nature of the dark energy. The candidates which have been proposed for explaining the 
dark energy are suitably described in a cosmological context which is generalized with respect to the 
ordinary one, admitting dark matter-energy couplings as well as generalized theories of gravity. 
In particular for the latter class of theories, a systematic treatment of the weak lensing process lacks in 
literature and this work aims at filling this gap. We considered a Lagrangian where the gravitational 
sector is made of a function which depends arbitrarily on the Ricci scalar as well as on a scalar field; the 
most general scalar-tensor theory of gravity, as well as the most general dependence on the Ricci scalar 
without a scalar field, can be described in full generality in this framework. 

We studied the generalized Poisson equations linking the fluctuating components to the two gauge 
invariant generalized gravitational potentials representing the metric fluctuations which cause the weak 
lensing process itself. This allowed to fix the contributions from the modified background expansion 
as well as the fluctuations, both in the Ricci scalar and in the scalar field responsible for the scalar- 
tensor coupling. We show that both of them are responsible for an anisotropic stress of scalar origin, 
causing the gravitational potentials to be different already at the linear level. We studied in particular 
the modifications induced by the time variation of the effective gravitational constant, which are most 
relevant in non-minimally coupled models in which the gravitational Lagrangian sector is a product of 
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a function depending on a scalar field and on the Ricci scalar; we focus in particular on the Extended 
Quintessence (EQ) scenarios, where the scalar field playing the role of the dark energy and responsible 
for cosmic acceleration today possesses a quadratic coupling with the Ricci scalar. 

Starting from the equation describing the geodesic deviation, we derived the generalized expressions for 
distortion tensor and projected lensing potential, tracking the effects due to the time variation of the 
effective gravitational constant and the contribution of the anisotropic stress; in particular, we show how 
the latter yields a correction proportional to the correlation between the gravitational potentials. 
Finally, we specialized our results to the description of two quantities which are most relevant for 
observations, i.e. the lensing convergence power spectrum as well as the correlation between the lensing 
potential and the Integrated Sachs- Wolfe (ISW) effect affecting the total intensity and polarization 
anisotropies in the cosmic microwave background radiation. 

By considering again the particular case of EQ cosmologies, we worked out approximate expressions 
for the corrections induced by the time variation of the effective gravitational constant on the lensing 
potential, lensing-lensing correlation angular power spectrum, convergence angular power spectrum as 
well as lensing-ISW correlation. We showed that the order of magnitude of these effects is of the order of 
8ttG^c/)q , where £ is the coupling constant and <j>o is the present value of the dark energy field. It may be 
noted how such correction may be relevant even if the underlying theory is close to general relativity, i.e. 
if the Jordan-Brans-Dicke parameter ojjbd = 1/327tG£ 2 0q is large, depending on the relative balance 
between £ and 4>q. 

Despite of these interesting indications in the particular case of EQ cosmologies, the formulas we devel- 
oped here have great generality, allowing a direct interpretation of the modern weak lensing observations 
in the context of cosmologies with generalized theories of gravity. 
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